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ABSTRACT 


The development of anisotropy in an initially isotropic spectrum 
is studied numerically for two-dimensional magneto hydrodynamic turbulence. 

The anisotropy develops due to the combined effects of an externally imposed 
dc magnetic field cuid viscous and resistive dissipation at high wave numbers. 
The effect is most pronounced at high mechaniceU. and magnetic Reynolds numbers. 
The anisotropy is greater at the higher wave numbers. 
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1. Introduct ion 

In the last several years, a systaMtic theory of magnetohydrodynamic 
(MHD) turhiilaice has developed, but it has dealt largely with the case in which 
no mean dc magnetic field is present. The presence of a mean dc Biagnetic field 
renders MHD turbulence inherently smisotropic , and one can no longer ted^e 

« 

advantage of the contraction of the statistical description that the isotropy 
^ provides. Here, we report a beginning on the problen of describing MHD tur- 
bulent phenomena in the presence of the anisotropy introduced by a mean dc 
magnetic field. We proceed numerically, using a two dimensional, spectral- 
method, incompressible MHD code which has been modified to include the pres«ice 
of a mean dc magnetic field. 

It would, of course, have been preferable to have used three-dimen- 
sional numerics. However, codes with adequate resolution to study three dimen- 
sional MHD turbulence at high mechanical and magnetic Reynolds nvmbers are reure, 
and groups with adequate computer resources to study MHD turbulence with than 
systematically are reu*er still. It has been possible to represent in two 
dimensions several of the features which we believe to be central to the three 
dimensional case. 

There are relatively few published results of laboratory measurements 
of MHD turbulence. Five importeuit papers are due to Robinson, Rusbridge, and 
Saunders {I 968 ), Rusbridge ( 1969 ), Robinson and Rusbridge (1971), Zweben, Menyuk 
and Taylor (1979), and Zweben and Taylor (I 98 I/. The first three of these 
report measuretients on the Culham-Harwell ZETA toroidal Z-pinch, and the last 
two, measurements on the UCLA Macrotor Tokaraak. The most important two features 
that emerged from the magnetic fluctuation measur^nents on both machines were: 
(1) the single-point frequency spectra were broad-bsuid, extending from a few 
hundred kiloherz down to the lower limit of the frequency resolution, with a 


notable absence of spikes or peaks; and (2) the vave>number spectra were 
peaked sharply in the direction perpendicular to the mean field, which mani- 
fested itself as a large ratio of parallel to perpendicular correlation l^ights 
(^10). Both of these behaviors will be seaa below to aaaerge from the evolution 
of initially isotropic spectra in two dimensions. 

The presence of a similar anisotropy is esseitlal, in one form or 
another, to a theoretical derivation of the Strauss equations (Strauss 1976, 
Montgomery 1982). The Strauss equation ere a reduced set of MED equations, 
intermediate between two and three dimensions, which have found wide applica- 
bility to tokaimk dynamics and may have considerably broader utility. It is 
unclear Iww to proceed in their derivation if the strong anisotropy is not 
present initially, and it wo\ild be reassuring to know that an initially isotro- 
pic spectrum would relax to an appropriately anisotropic one. 

The code used here has evolved from the spectral nnethod techniques 
given by Orszag (1971) and Patterson and Orszag (1971). Fyfe, Joyce, and 
Montgomery (l977a,b) and Matthaeus and Montgomery (198O, I98I) have used 
earlier versions of the present code to study forced dissipative MED turbulence 
selective decay processes, and the evolution of the sheet pinch, all in the 
presence of periodic boundeury conditions. Orszag and Tang (1979) have used 
a similar code to study small scale effects in two dimensions, and Pouquet 
(1978) has reported closure calculations for the same geometry. A general 
review of two dimensional turbulence has been given by Kraichnan and Montgomery 
(1980). None of the above work addresses the case of a finite mean dc magnetic 
field. 

The outline of this present paper is as follows. In § 2, the 
dynamical equations are described: they are only a slight modification of 

those used previously. The computational technique is briefly remarked upon 
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in § 3* Results from the computations are presented in § U, 5» and 6. TSie 
most important effect consistently observed is the development of strong 
anisotropy in the spectra from isotropic initied conditions. The anisotropy 
develops towards the kinds of spectra which seem to prevail in the reported 
measurements in ZETA and Macrotor. The development of the anisotropy dep«ads 
upon the combined effect of the external dc magnetic field strength and the 
size of the Reynolds numbers in an unexpected, but ultimately simple, way. 

A simple model of the effect is offered in § 7i where the results are briefly 
summarized and further directions for research are indicated. 
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2 . Dynamical Equatione 

The incompressible a dissipative MHD equations in tvo dimensions are 

used. The magnetic field consists of a constant mean field part 

plus a time-dependent zero^ean turbulent part B * (B , B , O) = 7x(e^&). The 

X y 

magnetic vector potential is e a(x, y, t) so that the Coulcaab gauge is employed. 

For all VEuriables, d/dz = 0. The velocity field y = (v , v , 0) » Vx(ejj<) is 

expressed in tezms of a stream function 4^ « 4^(x» y, t) and has zero mean. In 

what has become a standard set of dimensionless variables, the vorticity 

U) = o)e = 7xy, so that w * Similarly, the electric current J = Je„, with 

** Z z 

V^a = -j. 

Magnetic fields sire measured in terms of the initial root mean square 
turbulent field strength B. Velocities are measured in units of the Alfven 
speeu corresponding to B. The dimensionless viscosity and resistivitv v and ri 
are the reciprocals of mechanical and magnetic Reynolds numbers, respectively. 
The simplest form of the dynamical equations is, in the two-dimensional 
geometry, 

1^ + y*7w = B-VJ + vV^w + (1) 

and 

If * . (2) 

Without the B^ terms, eqs. (l) and (2) become those considered 
previously. In a recent derivation (Montgomery, 1982), the Strauss (1976) 
equations were re-derived using a perturbation expansion of the full set of 
incompressible MHD equations in three dimensions, in powers of B/B^. It was 
necessary to assume in the derivation that the time derivatives 3/3t remained 
of 0(1 ), or that no zeroth-order population of Alfven *aves was present. (The 
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linearized solutions for the three-dimensional case, as for eqs. (1) and (2), 
is just a superposition of Alfven waves with angular frequencies o)(k) * 
where k is the wavenumber.) An objective of the present computation is to 
see how a spectrum evolves which does contain an initially isotropic spectrum 
of Alfven waves. It will be shown that the spatial dependence of such a spectrum 
on the parallel spatial coordinate x becomes progressively relatively weaker 
with time. 
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3. Computational Technique 

The essence of the computational method Is that all physical fields 
are expanded in truncated Fourier series. The Fourier coefficients are 
stepped forward in time and are saved at predetermined time steps to provide 
a history of the dynamical evolution vf the field variables. To be explicit, 
Fourier representations of u) and a (for example) are (see, e.g. , Fyfe et al, 
1977a, b): 


co(x,t) = Ej^o)(k,t )exp(ik*x) 

(3) 

a(x,t) = Zj^a(k,t )exp(ik*x) , 

(1*) 


where k = (k , k ) and k and k are integers. Thus the dimension of the 
square box is chosen for convenience to be 2 tt. The Fourier coefficients 
retained lie in the range k . = 1 < |k| < k , where k is essentially 

limited by available computer time. Limitations on the Reynolds numbers are 
provided by the requirements that v and n be large enough so that tne Fourier 
coefficients for |kj > k^^ are suppressed. 

The first test of the code was to ran it with the dissipative terms 
removed. Conservation laws ceui then be tested ard comparisons can be made 
between the computed behavior of the Fourier coefficients and the predictions 
of the (unphysical) absolute equilibrium ensemble theory (ste, e.g., Kraichnan 
and Montgomery, I980). 

The second step was to compare dissipative results from two runs, both 
having twice the other's value of Bud thus having twice the spatial 

resolution. The third step was to run several cases with varying values for 
V and n at a particular value of to determine the influence of the Reynolds 
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numbers. Finally, the last step was to run a number of cases vhich differed 
only in the value of thereby investigating the effect that varying mean 
magnetic field strength has on the dynamics. 

Initial Fourier coefficients were chosen so that they would be non- 
zero only within a given annulus in k-space. The w(k,t) were initially non-zero 
only for <. k £ kg, and the a(k,t) were non-zero only for k^ £ k £ kj^. Within 
their respective annuli, the values of the u)(k,0) were chosen so that all 
|v(k,0)j = lw(k,0)l/k were equal, and the a(k,0) were chosen so that all 
|B(k,0)| = k|a(k,0)| were equal. The phases of the a(k,0) and u)(k,0) were 
assigned randomly. The specification of the initial Fourier coefficients was 
completed by giving values for the magnetic energy = Ej^|B(k,t )|^/2 and the 
kinetic energy = rj^|y(k,t ) |^/2 at t = 0. 

The total energy is E = Eg + E^, the ’’cross helicity” is 
P = Zj^y(k,t) ’B*(k,t )/2. The mean square vector potential is A = Ej^|a(k,t ) |^/2. 
E, P, and A are significant quantities in the theory of two dimensional MHD 
turbulence in the absence of a mean field; they are the only known non- 
dissipative invariants which umain invariant under tnuication of the Foiirier- 
expanded (v = 0 * n) version of eqs. (1) and (2) with B^ ■ 0. Tf # 0, E 
and P still have this status, but A does not. We may define R = Eg/E^ as the 
ratio of the energies. 

About a dozen different sets of initial Fourier coefficients have 
been rsed for various runs. Many of the sets exhibited similar behavior. 

The number of different sets of initied. conditions presented in this paper 
has been kept to a minimum, but for the runs discussed, there are others 
unreported for which the behavior was similar. Details of the runs explicitly 


discussed here are collected in Tables 1 and 2. 
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4. Non-dissipative Tests 

To test the spectred code, several cases were run with v ■ n ■ 0, 

and the results compared with absolute equilibrium ensemble theory. Such 

results are, of course, unphysical, and should be regarded as preliminary to 

the dissipative results presented in § 5* 

Equilibrium ensemble theory for the * 0 case was given by Fyfe 

and Montgomery (1976); the case 0 is recovered by simply deleting the 

third "rugged" constant of the motion A. The prediction is a simple equiparti- 
o o 

tion, <lv(k)p> = <|B(k)p> = const., independently of it. This is true for 

all realizable values of <E> and <P>. The ratio <R> = <Eg>/<E^> * 1, in sharp 

contrast to the B^ = 0 case. There is no crowding to the long wavelengths as 

k as there is when B = 0. 

max o 

Non-dissipative computations were carried out with B^ = 0 and B^ * 1 
for different sets of initial coefficients. Time averages were made of phase 
functions of the Fourier ^efficients and these were compared with ensemble 
averages. The two should be equal to the extent that the system is ergodic. 

The results of a single set of Fourier coefficients will be presented 
here, call it set A. For run Al, B =0, and for run A2, B * 1. For both 
runs, k = 16, E = l.C', the time step was (256) The total number of time 

steps was 12,800. For these conditions, k^ * k^ = 3 and k^ * kj^ = 5* Time 
averages were performed over the last 11,520 time steps. 

The comparison between the numerical results and the ensanble pre- 
dictions is given in Table 1 and in Figures 1 through 4. In Table 1, the 
behavior of E, P, A, and the evolving ratio = ^^^v shown as a 1‘unction 

of time. The invariarits E, P, and A for B * 0 and E and P for B » 1 are 

o o 

conserved to a few percent in all cases. Non-dissipative runs for initial 
conditions "B"(Eg = E^ = 1, k^ = 5, k^ «= 8, k| » 10, k^ « 13) will not be re- 
ported in detail . 





4 





11 


In Figs. 1 through 4, the directionally averaged magnetic and kinetic 
modal energies are presented for time averaged data, (Directional averaging 
means averaging over all values of k corresponding to a particular k . ) At 
the lower values of k, where the number of degenerate Jt's is sparse, occasion- 
ally large anisotropies necessarily appear, but no systanatic directionality 
vas observed. An explicit comparison of these runs with corresponding dissi- 
pative runs will be presented in the next section, after suitable measures 
of anisotropy are defined. It will be seen that the non-dissipative runs, 
in contrast to the dissipative ones, show no anisotropy when averaged over 
long times. 

In Figs, 1 through the spectral predictions of the absolute 
equilibrium ensemble theory are shown as solid lines and the plotted points 
are the time-averaged results of the computations. No significant deparcures 
from the predictions of the absolute equilibriimi ensemble theory have been 
observed in these runs, and in others not reported here. Typiced. behavior 
for Imost all the individual Fourier modes (those with k #0) is that of 

X 

an Alfven wave of angular frequency to(k) = ik*B . This is the principal 

- -- -»o 

qualitative difference in the time behavior observed between the = 0 and 
B^ 5 ^ 0 cases. Superposed on the Alfven-wave oscillations is a slower transfer 
of excitation between the Fourier modes. 

In the next few sections, we pass to a consideration of the dissipa- 
tive cases. In order to maintain a common tliread throughout the discussion, 
we will primarily use set A as the initial conditions for the dissipative 
runs to be presented. The principal quantities to be varied are the spatial 
resolution the values of the dissii>ation coefficients v and n, and 

the external field strength B . 
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5 . Appearance of Anisotropy; Effects of Viscosity, Realstlvlty, and Spatial 
Resolution 

In this section, ve describe dissipative rvns. The pres^ce of 

dissipation is always central in the evolution of any real turbulsnt field. 

We use primarily the saroe initial Fourier coefficients as those in 

§ U. The most striking effect ve have observed is that when both dissl|«tion 

and a mean B a'.e present, anisotropy appears. The k spectrum e/olves into 
o 

one peaked perpendicularly to the mean field F . We describe the dependence 
of this anisotropy on the dissipation coefficients euid the mean B^, after first 
r«narking upon some considerations of necessary spatial resolution. 

For an accurate solution to eqs. (l) and (2), we must resolve the 
smallest dynamically significant spatial scales that v and n permit. An estimate 
for these is provided by the Kolmogoroff ’’dissipation wave number”, constructed 
by dimensional analysis based on the rate of dissipation of energy. For MHD 
in two dimensions, it is 

k^ - [n“^|dydt.|^-^ v“3|dE^/dtj^]^''\ (5) 

where [dEg/dtl^^ and |dE^/dt|^ are the ohmic and viscous energy dissipation 
rates, respectively. Ideally, k^ should be less than for accurate solu- 
tion of eqs. (1) and (2) eni this provides at present the most severe limita- 
tion on turbulence computations, both Navier -Stokes and MHD. For son» 
purposes, it may be that when the smallest scales ai’e dynamically insignificant, 
kj somewhat lerger tlmn k may be tolerated, but k.»k must always I'ignal 

u XmX Cl 9lu 

uselessness In a turbulence computation. Other measures of m<^ Iwgth scales 
for the turbulent field, such as those of Pao and T^lor (see, e.g. , Leslie 
19T3) may be used to characterize turbulent activity, but the iCslmogoroff 
scale is the one most widely accepted and ve shall use it hereafter. 



The dissipation wave number may be adjusted empiricaLlly, for a given 
set of initial Fourier coefficients, by raising v and n until the computed 
remains no greater than . We consider here only the case v « n (unit 
magnetic Prandtl number). 

In the initial set of dissipative cases, there are six runs (A3 
tlxrough A8, see Table 2 ), having = 0 and = 1 for v = ri = 0.005, 0.01 
and 0.02. The following characteristics were conroon to these six runs: 

k =32, time step = (256)”^, total number of time steps = 1280, initial 

max 

Fourier coefficients: set A. The tonporal evolution of the variables E, P, 

A, and for the runs with v = n = 0.01 are shown in figures 5 and 6. (Qv^i- 

tatively similar behavior was observed for v = n = 0.005 and 0.02, the decay 

times simply increased and decreased, respectively.) 

When = 0, roughly twice as much magnetic as kinetic energy 

develops; when = 1, the ratio remains approximately unity. E, P, A all 

monotonically decay for B = 0, but when B = 1, A actually increases at times, 

o o 

but overall decays. 

To measure anisotropy in the k spectrum, a set of angles were defined 
for each Fourier-decomposable field, 

= tan-^ { = (6) 

where Q is any one of the fields y, w, a, B, or J. For a k spectrum purely 

normeil to B^, 0^ is 90°, and for an isotropic spectrum, 1*5°. 

Temporal evolution of the angles for B = 0 emd B * 1 is illustrated, 

o o 

for runs A^> and A6, in Figs. 7 and 8, which both refer to V = n = 0.01. For 
= 0, the angles roam unsystematically about ^5*^, indicating isotropic 
developnent . The anisotropy in Fig. 8, for = 1, is typical. Notice that 


Ik 


the angles tend to increase in the following order: ^ ^ 

6^ < 6g < 0^: the anisotropy is most prono^mced at the shorter wavelengths. 

Also, when we consider v = n * 0.005 and 0.(S as well, we see tl»t the degree 
of anisotropy tends to decrease with increasing V and n> The results are 
summarized in Table 3, which displays averages of the 6 q for v ■ n * 0.005* 
0.01, and 0.02 between times 1.5 and 5*0. 

Also calculated were the dissipation wave number (eq. (6)), mean 

2 2 
square vorticity = Ij^|ui(k.t ) | , and mean square current J = I^|j(k,t)Pi 

the lott-er two quantities are required to evaluate eq. (6). The taaporal 

evolution of Q and J for runs A5 and A6 (v * n = 0.01) is shown in Fig. 9. 

Table It displays the maximum values of k^, fl, and J for the six runs A3 

through A8. In Table It, it is apparent that the maximxan dissipation wave 

number decreases with Increasing v and n » as do the maxima of J and fi. The 

presence of a finite also apparently enforces a more nearly equal partition 

of dissipation between 0 , and J. 

From the foregoing results, it is clear that the maximum k^ which 

will fit inside the resolution k =32 occurs for runs A5» A6, where 

max 

V = n = 0.01. To test the effect of lowering the resolution below k on 

max 

the results of runs A5 and Ab, runs A9 and AlO were made with identical para- 
meters to A5 and a 6, tut with k = l6. (¥e also showed that halving k sd- 

max max 

lowed doubling the time step. ) In runs A9 and AlO, the time histories of the 

angles 0^ and the quantities k^» fl, and J were essentiedly the same as those 

for A5 and A6 (see Shebalin, 1982, for further details). This insensivity 

to halving the spatial resolution, so that k ^^^ falls well inside k^, gives us 

seme confidence in the validity of the results of such runs at the higher 

Reynolds numbers, such as runs A3 and AU, where k lies outside k, . 

IP ft X ^ 
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We should remark at this point on the question of the anisotropy of 
the non -dissipative runs. In Fig. 10. ve display the evolution of the angles 
0^ and 0g for a dissipative run B2 (v * n * 0.0025). Both nms have the ssjm 
initial Fourier coefficients with = 5» ® 8» k^ * 10. k^ * 13. Eg “ for 

both rvins, and the non-vanishing Fourier coefficients have equal aapliti«ies 
within their respective annuli, and random phases. (See Table 2.) 

Referring to Fig. 10. we can see that for Bl. although 0„ and 0_ 
initially rise to about 51° and 55° respect iv^y. they subsequently fall and 
oscillate about ^5°. Run B2 , however, shows 0^ and 0g rising to maxima of 
approximately 66° and 68°. respectively, and then oscillating near those values. 
At the beginning of the dissipationless run, there apparently is a t«idency 
toward anisotropy which cannot maintain itself; the v = n = 0.0025 run raaains 
strongly anisotropic, however. This behavior, involving the necessity of 
small but finite dissipation in the maintenance of anisotropy, corresponds, 
we believe, to a simple physical effect which is discussed in § ?. It is our 
belief that no long-time anisotropy is to be expected without the presence of 
dissipation, and (for reasons given in § 7 ) that the smaller the dissipation 
coefficients, the greater the degree of anisotropy is likely to be. A computa- 
tion which could afford the spatied resolution to resolve signifiwantly 
higher values of k^j than we can resolve would see proportionately higher 
degrees of anisotropy, we believe. 


16 


6. Effects of Variable Mean Field Strgagth 

Having discussed the effects of the spatial resolution and 

the dissipation coefficients V and n on the results, we pass to a consideration 

of the effects of varying the mean field Set A is chosen for the initial 

Fourier coefficients, and k = l6. Thus the runs in this section are initially 

max 

similar to A9 and AlO, except that takes on the values l/l6, 1/8, l/^t, 1/2, 

2, 4. 8, and 16. These values correspond to runs All through Al8, respectively. 
(Parameters for all runs considered appear in Table 2.) 

Rather than display a multitude of graphs similar to those presented 
for rvins A9 and AlO, we display mainly graphs of time averaged angles 0 q. In 
Fig. 11, we show 0^, averaged over times 1.5 to 5-0, as functions of B . In 

(qj O 

Fig. 12, we show, as fxinctions of B^, the maximum value of k^, the time when 
this maximum occurs, the maximum enstrophy and the maximum mean square 
ciirrent J. 

Figures 11 and 12 illustrate a number of interesting effects. First, 

there is anisotropy which develops as increases from zero. Ser-ond, the 

effect saturates: beyond a value of B of about 2, further increase in B 

o o 

results in no further increase in anisotropy, for these values of k , v, and 

max 

n. Third, the values of 71 and J approach each other as increases, reflect- 
ing a progressively more Alfven -wave-1 ike behavior at the dissipation scales. 

The anisotropy continued to be most pronoxuiced at the highest wave niimbers: 

6 Eind 0 were larger than the other angles. 

'A) j 

A physical feeling for the configuration-space manifestations of 

the anisotropy may be obtained from Fig. 13. There, contour plots are given 

for the vort’ -ity w(x,t) and current J(x,t) at various times with a zero emd 

a non-zero value of B . For the plot in which B = 2 it is clear that, at 

o o 
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time t = 2.0, the vortlcity and current co.''tours have elongated in the direc- 
tion of the mean field, raainiscent of the elongations that vere reported in 
the ZETA and Macrotor devices. 


» 


c 


r 
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T. Discussion 

The linearized MHD equations in tvo and three dimensions result 
from discarding the right hand sides of 


- B -7B - vV^v = B-VB - v7v - Vp (7) 

3t~o 

and 


2 

— - B *Vv - nv B = B-7v - vVb, (8) 

dt ~0 ~ ■“ - ~ - 


subject to V’v = 0 = 7*B. Temporarily ignoring the dissipation, the most 
general linearized solutions can be vritten as B = b^^ + y - -b^ + bj^. 


where 


b = E bfexp(ik*x - io)(k)t) + c.c. 


(9) 


b = E b exp(ik*x + ioj(k)t) + c.c. 

•vjj K K -w 


with o)(k) = "c.c.'* standc for complex conjugate, and the amplitudes 

R L R L 

b * satisfy k*b^’ = 0, but are otherwise arbitrary. The division of the 
fields, of course, is into right and left traveling waves. 

Assuming that eqs. (9) are a satisfactory zeroth order solution to 
eqs. (7) and (8), we may inquire, within a perturbation-theoretic framework, 
about the effects of directionality on the modal transfer. We may substitute 
the linear solutions into the right hand sides of eqs. (7) and (8) and proceed 
iteratively to calculate the first nonlinear correction to the linear fields. 
Detailed inspection shows that there is no net coupling between the right- 
traveling waves with each other or the left traveling waves with each other. 
The only non-zero couplings are between right and left traveling waves. 


19 


To resonate effectively with a third, initially-unexcited Fourior 
mode, there sure matching conditions on both frequaicy suid wavenumber which 
must be met. These are very restrictive whai it is taJEen into consideratitHi 
that one interacting wave must be right traveling and the other left traveling. 
If the two waves have wavenumbers ^ and frequencies - (*>(kg), the 

conditions that they be able to excite a third wave resonantly with wavenumber 
k^ and frequency ±u>(k 2 ) are that 

±u)(k2) = w(k^) - oiCkg). (10) 

Since u)(k) = eqs. (10) have a solution only if either w(^) = 0 or 

w(k 2 ) = 0, so that either cr ^2 has zero component along B^. Thus, a 

three-wave resonsint interaction can result in the excitation of a wave with a 
larger value of [k^l than that of either of the other two, but never with a 
larger value of lk^|. It is clear that excitations may readily transfer 
energy by this process in the perpendicular direction in k space but not in 
the parallel one. 

An initisJ-ly isotropic distribution in k space elongates in the 
perpendicular direction vuitil something stops the migration to larger [k |. 

y 

In the present computation, that is either k^ or, for the truncated non- 

dissipative model, k . In the latter case, eventusd isotropization occurs 

max 

as a consequence of higher order processes. In the presence of finite dissipa- 
tive decay, the anisotropy persists, as in Fig. 10. 

We believe the above-described mechsuiism to be responsible for 
the observed anisotropy of magnetic fluctuations in toroidal devices. For a 
given level of excitations, k^ increases as v and n decrease, so the effect 
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should he most pronounced at high Reynolds numbers. The independence of 
further increases in beyond a certain nK^dest level, must simply mean that 
transfer in the parallel k-space direction has been effectively frozen relative 
to the transfer perpendicular to B beyond a value of B^ Just a few times the 
mean fluctuating field strength. 

The natural directions into which these investigations shovild be 

taken are to some degree obvious. First, considerably higher Reynolds numbers, 

with their necessarily higher spatial resolution, should be investigated. 

Second, spatially inhomogeneous vacuum fields shovild be added to ascertain 

the effects of mean field curvature on the anisotropy. Finally, the effects 

of the mean fields on the small scales need to be ascertained: how does B 

~o 

affect current filamentation, x-point behavior, and magnetic re-connection? 
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Table 1 . Non-dissipatlve runs (v « 0 ■ ri) tor comparison vith absolute 
equilibrium ^isemble theory [Tia« step: (2 56}**^. Total no. of time steps: 
12,800. Averages taken over last 11»520 time steps.] 



Initial 

Final 

Tiro 

Average 

% Change 
(Initial -Final) 



A1(B - 0) 
o 



E 

1.000 

0.97li* 

0.9782 

2.9 

P 

0.235*^ 

0.2167 

0.2227 

7.9 

A 

0.03288 

0.03286 

0.03288 

0.06 



A2(B^ « 1) 



E 

1.0000 

1.0150 

1.0062 

1.5 

P 

0.235^ 

0.2277 

0.2313 

3.3 

A 

0.03288 

O.OIU97 

0.01489 

54.5 


Initial 

R - Eg/E^ 
Final 

Time 

Average 

Theoretical 

A1 

1.0000 

1.0058 

1.0608 

1.0484 

A^ 

1.0000 

0.9825 

1.0004 

1.0000 


Table 2. Run Parameters 

If 

Run max 

ORIGINAL PAGE 13 
OF POOR QUALITY 

v*n ®o 

Time Step Size 

2k 

Total 0 of 
Time Steps 

A1 

16 

0 

0 

1/256 

12600 

A2 

16 

0 

1 

1/256 

12800 

A3 

32 

.005 

0 

1/256 

1260 

AU 

32 

.005 

1 

1/256 

1280 

A5 

32 

.01 

0 

1/256 

1280 

A6 

32 

.01 

1 

1/256 

1280 

AT 

32 

.02 

0 

1/256 

1280 

A8 

32 

.02 

1 

1/256 

1280 

A9 

16 

.01 

0 

1/128 

61*0 

AlO 

16 

.01 

1 

1/128 

6 U 0 

All 

l6 

.01 

l/l6 

1/128 

6k0 

A12 

16 

.01 

1/8 

1/128 

6k0 

A13 

16 

.01 

1/i* 

1/128 

6k0 

All^ 

16 

.01 

1/2 

1/128 

61*0 

A15 

16 

.01 

2 

1/256 

1280 

A16 

16 

.01 

k 

1/512 

2560 

A17 

16 

.01 

8 

I/I 02 L 

5120 

A13 

16 

.01 

16 

1/201*8 

IO 2 U 0 

B1 

16 

0 

1 

1/128 

3200 

B2 

16 

.0025 

1 

1/128 

3200 

Set A: 

initial equipartion 

of ener^ 

in all magnetic and kinetic Fourier modes. 

Set B: 

initial equipartion of energy 

in Foxu*ier modes such that 9 

1 1 25; 

Initi6Ll 

Fourier modes outside this annulus initially set to zero, 
ms values of and |v|^ were both equal to one for sets 

A and B. 
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Table 3» Average Angles between t ■ 1.5 and 5-0 


v.n 

A. 

e 

a 

e 

V 

A 

6 

h) 


.005 

53.4 

50.6 

58.2 

59.3 

64.0 

67.5 

.01 

51.6 

48.5 

54.8 

55.0 

60.5 

64.0 

.02 

48.7 

46.5 

50.8 

50.2 

53.6 

56.6 
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Table U . 

Maximum Values 

of K 

( 

Run 

v^n 

B 

0 

A3 

.005 

0 

AU 

.005 

1 

A5 

.01 

0 

a6 

.01 

1 

A7 

.02 

0 

a8 

.02 

1 


fl and J for runs A3 through A8 


Max 

Max 

Max 


Q 

_J 

49.6 

25.0 

51.8 

48.3 

30.6 

38.0 

32.2 

IT.O 

39.0 

31.4 

21.4 

27.8 

21.0 

17.0 

26.3 

20.8 

17.0 

20.3 
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fimm cAPTioKg 


Fig. 1 . Direct everegtA klnttie wergy epectre tor rm Ali- 

B ■Ot\^*n*0tk ■ Solid line it theoretical j^edietic». 

o ttax 

_2. Directionally averaged aagnrtie spectra for run Al! 

■ 0^ V • n ■ 0» « 16. Solid line Is theoretical prediction. 

J[, Directionally averaged kinetic energy spectra for rim Ml 

B ■ I, V ■ n * 0» k ■ 16. Solid line is theoretical predict!^, 

o hiax 

Directionally averaged wigretie energy ipeetra for rm Mi 

B ■ I, V ■ n ■ 0, k ■ 16, Solid line is theoretical {prediction, 

o wax 

The quantities E, P» A* and R for run A5; ■ 0, v»ii ■ O.Ol* 

k • 3£. 
max 

The quantii.'^'i E, P, A» and R for run A6: ■ 1, v ■ H • 0.01, 


Fig. 4 . 


Fig. 5 . 


Fig, 


EMr. .1* 
Fig, 


Fig. 10 . 

Fig. 11. 
Fig, la . 


wax 

Measures of anisotropy for run A5; ■ 0, v ■ n ■ 0.^» " 32. 

Measures of anisotropy for run A6: ■ 1, v ■ n ■ 0.01, It . • 3^. 

a) Dissipation wave mmber for run A5* b) dissipation wave nimber 
for run A6j c) mean square vortlclty for run A5» d) man square cur- 
rent for run A5i e) wean square vortleity for run A6; f) wean aquare 
current for run Ab. 

The angles 6 and 0„ for runs B1 (B^ ■ 1, u ■ p ■ 0) and B2 

V B O 

(B » 1, V • n ■ o.oogj). 

o 

Measures of anisotropy at a function ©f mean field strength B^. 

a) Maxlrnw k. at a function of B t b) time of wax k.; c) mxlww me«n 
a o a 

square vortlclty at a function of d) mxlmim mean square current 

as a function of B . 

o 
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Fig. 13 « a) Vorticity of set A at t = 0; b) current of set A at t « 0; 

c) vorticity of run A9 (B^ « 0, v ® n ® .01, ® l6) at t = 2.0; 

d) current of run A9 at t = 0; e) vorticity of run A15 (B = 2, 

o 

V = ri = .01, k = l6) at t = 2.0; f) current of run A15 at t = 2.0. 
max 

t 


( 


♦ 


IT 
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16. Solid line is theoretical prediction. 
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l6. Solid line is theoretical prediction. 


OF POOR QUALIYY 


* 


f 





16. Solid line is theoretical prediction. 


npieiNM. 


o 

u*^ 



o 

CD 

UD 



O-J 

CD 

PmI 



CD 


o 

<s» 

CD 


iS 

O 

(S*l 

iH 

<St 

C9 

O 



CD 

SI 


o 

CD 

CD 





tcQ 


R 



* 

O 

R 

CT 

R 

> 


H 

R 


w 


l6. Solid line is theoretical prediction. 
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Fig. 3 . The quantities E, P, A, and R for run A5: = 0, v=n * 0.01 
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Fig. 6 . The quantities E, P, A, and R for 
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Fig. 7 » Measures of anisotropy for run A5: B = 0, v = n = O.ttL 



Measures of anisotropy for run A6; B = 1, v * n * 0.01 
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